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The thermoelectric properties of two-dimensional graphyne sheets are investigated by using first- 


principles calculations and the Boltzmann transport equation method. The electronic structure indicates 


a semiconducting phase for graphyne, compared with the metallic phase of graphene. Consequently, 


the obtained Seebeck coefficient and the power factor of graphyne are much higher than those of 


graphene. The calculated phonon mean free path for graphene is 866 nm, which is in good agreement 


with the experimental value of 775 nm. Meanwhile the phonon mean free path of graphyne is only 


60 nm, leading to two order lower thermal conductivity than graphene. We show that the low thermal 


conductivity of graphyne is due to its mixed sp/sp* bonding. Our calculations show that the optimized 


ZT values of graphyne sheets can reach 5.3 at intermediate temperature by appropriate doping. 


1 Introduction 


As thermoelectric materials can directly convert heat to elec- 
tricity and vice versa,’ they have been attracting increasing 
attention from both theoretical and technical sides. The conversion 
efficiency is governed by the dimensionless figure of merit 
ZT = S’oT/k, where S, c, and T are, respectively, the Seebeck 
coefficient, the electrical conductivity, and the absolute temperature. 
The thermal conductivity x includes electronic (xe) and phonon 
contributions (xpn).” To get high ZT values is usually difficult 
since the transport coefficients S, o, and « are coupled with each 
other. However, it was predicted that these factors could be 
decoupled in low dimensional systems.’ As a representative of 
low-dimensional carbon materials, graphene has drawn wide- 
spread interest in many fields. The electron mobility of graphene 
is up to 15000 cm* V's‘, which is a great advantage to 
thermoelectric applications." Unfortunately, graphene has a 
very high thermal conductivity of 2500-5000 W m~t K*,°° 
leading to a low ZT value of 0.02.° 

Graphyne is a new type of two-dimensional carbon allotrope. 
Unlike the sp* bonded graphene, graphyne contains two kinds 
of bonds, namely sp and sp”. Recently, large graphyne film has 
been successfully prepared and it exhibits excellent semicon- 
ducting properties.*° It is conceivable that graphyne may be 
more suitable for thermoelectric applications than metallic 
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graphene. Indeed, more and more attention has been paid to 
the electronic, thermal and thermoelectric properties of graphyne. 
With respect to its electronic properties, density functional 
theory calculations indicated that a band gap is opened in 
graphyne by the presence of the sp bonding.**’* Shuai et al. 
calculated the intrinsic mobility in graphyne using first- 
principles calculations, and it was found that the carrier mobility 
in graphyne is even larger than that in graphene.’*’” In the 
aspect of thermal conductivity, Zhang et al.‘* have shown, based 
on their molecular dynamics (MD) simulations, that the thermal 
conductivity of graphyne is one order of magnitude lower than 
1.°° applied the 
nonequilibrium Green’s function (NEGF) approach to graphyne, 
and it was found that the thermal conductance of graphyne is 
40% of that of graphene. By using a similar method, Ouyang 
et al.*°”! have shown that the thermal conductance of B- and 
y-graphyne nanoribbons is 26% and 40% of that of graphene 


that of graphene at room temperature. Lu et a 


nanoribbons, respectively. In the aspect of thermoelectric properties, 
Wang,” Ouyang,” and Sevinçli” predicted that the optimized 
ZT values of graphyne sheets can be one order higher than that 
of graphene, mainly owing to the increase of the power factors. 
Sun et al.’ carried out MD simulations and showed that the 
lattice thermal conductivity of graphdiyne is 7.3 W m™* K * at 
room temperature, resulting in a ZT value of 4.8. 

In previous studies, the NEGF method at the ballistic limit is 
employed to investigate the electron and phonon transport 
properties. It is known that electron-phonon and phonon- 
phonon interactions are not well treated in the NEGF method.”* 
Some important information, such as the electron or phonon 


mean free path, can neither be extracted from the calculations. 
The Boltzmann transport equation (BTE) coupled with deformation 
potential theory (DPT) was applied to calculate the carrier 
mobility of some carbon materials t??? and this method 
was proved to be suitable for describing the electron-phonon 
coupling of graphyne. Meanwhile, an exact numerical solution 
of the phonon BTE has been developed to study the thermal 
conductivity. This method has been successfully applied in the 
phonon transport study for two-dimensional materials such as 
graphene®”® and WS3.” In this work, we use the BTE method 
based on first-principles calculations to investigate the electronic, 
phonon and thermoelectric properties of graphyne. By considering 
the phonon-phonon scattering, our theoretical calculations 
find that the thermal conductivity of graphyne is two orders 
lower than that of graphene, leading to the much higher 
thermoelectric performance than the previous predictions. 

The rest of this paper is organized as follows. Section 2 gives 
the computational details of our theoretical methods. In Section 3, 
we discuss the electronic, phonon, and thermoelectric transport 
properties of y-graphyne and graphdiyne. A summary of our 
work is given in Section 4. 


2 Computational details 


The calculations of graphyne are performed using a plane-wave 
pseudopotential formulation???” within the framework of density 
functional theory (DFT). The generalized gradient approximation 
(GGA) of Perdew-Burke-Ernzerhof (PBE)? is adopted for 
exchange-correlation energy. Projector-augmented wave (PAW) 
potentials are used for the carbon atoms and the cutoff energy is 
set to 400 eV. The Brillouin zone is sampled with 20 x 20 x 1,14 x 
14 x 1, and 12 x 12 x 1 Monkhorst-Pack &meshes including the 
T-point for graphene, y-graphyne and graphdiyne, respectively. A 
vacuum layer of 10 A is added along the z axis for the two- 
dimensional structures. During the structure optimizations, 
both the atom positions and lattice constants are fully relaxed 
until the magnitude of the force acting on the atoms is less than 
0.001 eV A~*, which also converges the total energy within 1 meV. 

The electronic transport coefficients are obtained from the 
analysis of band structure calculations by using the BTE** 
method within the rigid band approximation.’ This semi- 
classical approach has been successfully applied to a wide 
range of thermoelectric materials.*° The kernel is to find the 
conductivity tensor 
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where N is the number of k-points sampled. To get reliable 
results, a very dense k mesh up to 4000 points for graphyne is 
used in BTE calculations. The k-dependent conductivity tensor 
Gaplik) is given as 


Oup(ik) = eT iWValik)vglik), (2) 


where q(ik) is the carrier relaxation time and the group velocity 
is v(ik) = (O&,/0k)/i. The electrical conductivity o, the Seebeck 


coefficient S, and the electronic thermal conductance x, can be 
calculated by integrating the conductivity tensor** 
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where u is the chemical potential, f,(T,¢) is the equilibrium 
Fermi-Dirac distribution, and Q is the volume of the unit cell. 
Usually, Q is defined for three-dimensional systems. As to two- 
dimensional graphene and graphyne sheets, a layer thickness is 
required to define the electrical and thermal conductivity. In 
order to connect to the previous studies, this value is conven- 
tionally chosen to be ô = 3.4 Å, which is the separation between 
graphene layers in graphite. 

The carrier relaxation time t in eqn (2) measures how 
quickly the carrier can restore the equilibrium state after 
scattered by phonons. The DPT method proposed by Bardeen 
and Shockley?” can be used to describe the electron-phonon 
scattering in the long-wave limit. In two-dimensional systems 
like graphyne, the relaxation time is given byt’ 
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where the parameters ñ and kg are the Planck constant and the 
Boltzmann constant, respectively. In addition, m* is the density 
of state effective mass, C is the elastic constant, and Epp is the 
deformation potential constant. These three parameters can be 
readily obtained from first-principles calculations. Note here m* 
contains contributions from ky and k, directions, and the degeneracy 
of light and heavy bands needs to be taken into account. 

We calculate the phonon thermal conductivity of graphene, 
y-graphyne and graphdiyne by solving the linearized BTE 
numerically as implemented in the ShengBTE code.” A temperature 
gradient across a material results in a heat flow because of the 
phonon diffusion from hot to cold. However, these phonons 
undergo scattering in the diffusion process and the phonon 
distribution f,, is changed, where s is the branch index and q is 
the wave vector. The phonon distribution is affected by two 
factors: diffusion from the temperature gradient and scattering 
due to the allowed processes. In the steady state, the rate of 
change in the phonon distribution is vanished and it can be 
expressed by the BTE:*? 
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specific scattering processes, including phonon-phonon scatter- 
ing, impurity scattering, defect scattering, and so on. Practically, 
when the only scattering sources are the two- and three-phonon 
scattering processes, a linear version of the BTE is adopted and 


solved for scattering time tqs by using an iterative method.*? 
Finally, the phonon thermal conductivity is given by the formula:° 


T agal (hog) f gs (Sas + 1)¥qs qed. (8) 
Here, @gs and vgs = Õwgs/ðq are the phonon frequency and the 
group velocity, respectively. 6 = 3.4 A is the layer thickness of 
graphene and graphyne sheets as mentioned before. 

Based on the second- and third-order interatomic force 
constants (IFCs) calculated using first-principles calculations, 
the phonon dispersions and scattering properties are deter- 
mined. The results of phonon thermal conductivity are care- 
fully tested with respect to the size of the supercell, the nearest 
neighboring range, and the q-point mesh. In the calculation of 
second- and third-order IFCs, we use 10 x 10 x 1and3 x 3x 1 
supercells for graphene and graphyne sheets, respectively. The 
interactions up to the fourth nearest neighbors are considered 
when dealing with the anharmonic one. The g-point mesh is set 
to 100 x 100 x 1, 40 x 40 x 1 and 30 x 30 x 1 for graphene, 
y-graphyne and graphdiyne, respectively. 


3 Results and discussion 
3.1 Crystal and band structures 


The structures of two-dimensional y-graphyne and graphdiyne 
are illustrated in Fig. 1. These two kinds of graphyne sheets 
could be considered as benzene rings connected by carbon 
triplet bonds. The structural symmetry of y-graphyne and 
graphdiyne is P6/mmm, the same as graphene. As shown in 


(a) HALOTT 


~N 
© 
o 


Fig. 1(b), graphdiyne contains one more acetylenic linkage between 
the nearest neighboring carbon hexagons than y-graphyne 
(Fig. 1(a)). The optimized lattice constants ao are 6.890 and 
9.460 Å. The bond lengths depicted in Fig. 1 are as follows: for 
y-graphyne, 1, = 1.426 A (sp?-sp*), J, = 1.408 A (sp-sp°), and 
l = 1.223 A (sp-sp), and for graphdiyne, L = 1.432 A (sp?-sp”), 
l, = 1.396 A (sp-sp’), and J; = 1.233 A (sp-sp), and l, = 1.339 A 
(sp-sp). The obtained structural parameters are in good agree- 
ment with previous LDA** and GGA calculations.’”"**° 

Fig. 2(a-c) presents the calculated energy band structures of 
graphene, y-graphyne and graphdiyne. As is known, the electronic 
structure of graphene is characterized by the Dirac cone, where 
the conduction and valence bands meet at the Dirac point. In 
y-graphyne and graphdiyne, the band crossing at the Dirac 
point is broken due to the presence of sp bonding, showing a 
semiconducting behavior. As shown in Fig. 2b and c, y-graphyne 
has a direct band gap of 0.46 at the M point, and graphdiyne has 
a direct band gap of 0.48 eV at the I’ point. The conduction band 
maximum (CBM) and valence band minimum (VBM) of y-graphyne 
are singly degenerate, while those of graphdiyne are both doubly 
degenerate. The band structures presented in Fig. 2 are consistent 
with the previous studies.” 

Using the formula m* = fi’/(0°E/0k), the effective masses at 
CBM and VBM can be calculated. As listed in Table 1, the 
effective masses in y-graphyne are anisotropic - the effective 
masses along the F —> M direction are much larger than those 
along the K > M direction. In contrast, the effective masses in 
graphdiyne are rather isotropic. For both y-graphyne and 
graphdiyne, there is no significant difference between the 
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Fig. 2 Calculated band structures of (a) graphene, (b) y-graphyne, and (c) graphdiyne. The Fermi level is at O eV. 


Table 1 Band gap E, and effective mass in the conduction band (m,*) and 
the valence band (m,*) for y-graphyne and graphdiyne 


Structure Eg (eV) m.*/mMo m,*/Mo 

y-Graphyne 0.46 at M 0.197 0.082 0.215 0.085 

Ref. 11 0.52 atM 0.15 0.063 0.17 0.066 

Ref. 43 0.46 at M 0.20 0.080 0.21 0.083 

Graphdiyne 0.48 at I 0.088 0.090 0.089 0.097 
(Mor) C-K) (MoT) (eR 

Ref. 11 0.53 at I 0.073 0.075 


effective masses of holes and electrons. Narita et al." and Kang 
et al.” investigated the band structures and effective masses of 
y-graphyne and graphdiyne using local density approximation 
(LDA) and GGA methods. Our calculated results are consistent 
with their results, as summarized in Table 1. Except for the 
slight difference in numerical values, GGA and LDA present 
similar band gaps and similar effective masses. 


3.2 Electronic transport coefficient 


Based on the calculated band energies in the first Brillouin 
zone, the electronic transport coefficients of graphyne can be 
evaluated by using the semi-empirical BTE approach with 
the rigid band approximation. Within the BTE method, the 
electrical conductivity ø, the electronic thermal conductivity Ke 
and the power factor So are calculated with respect to q, that is 
to say, the obtained electronic transport coefficients are o/t, Ke/t 
and S’o/t, respectively. The realistic relaxation time t is energy 
dependent. However, it is often reasonable to assume a constant 
t since the electrical conductivity determined by t does not 
change very much within the energy scale of kgT. In this work, 
the relaxation time is calculated by applying the DPT method 
with the effective mass approximation. 

Besides carrier effective mass, the two other parameters in 
eqn (6), the elastic constant C and the deformation potential 
constant Epp, can also be readily obtained from first-principles 
calculations. Specifically, we change the lattice constants of 
graphyne by +£0.5%, +1.0% and +1.5% along the a and b axes, 
and then measure the variation of the total energy and the shift 
of the VBM and the CBM. The two-dimensional elastic constant 
C is obtained by fitting the total energy E with respect to the 
dilation Al/Ip, as (E — Eo)/So = (C/2)(Al/I)”, where So is the area 
of the cell in the plane and Eo is the total energy in the 
equilibrium state. The deformation potential constants Epp 
for electrons and holes are obtained by fitting the shift CBM 
and VBM AE with respect to the dilation, as Epp = AE/(Al/lọ). 
Using the obtained effective masses, elastic constants and 
deformation potential constants, the relaxation times can be 
calculated according to eqn (6). 

The calculated relaxation times at room temperature as well 
as two-dimensional elastic constants and deformation potential 
constants for y-graphyne and graphdiyne are summarized in 
Table 2. Shuai’s group investigated the carrier relaxation time in 
carbon nanomaterials including graphyne.'*'” Our calculated 
results are very much consistent with their calculations, as 


Table 2 Elastic constants C, deformation potential constants Epp, and 
carrier relaxation times t for y-graphyne and graphdiyne at 300 K 


C EBp Epp q" g 

Structure Axes (Jm’) (eV) (eV) (ps) (ps) 
y-Graphyne Present work a 220.35 4.56 2.81 0.49 1.41 
b 220.38 4.57 2.79 0.49 1.43 
Ref. 17 0.56 1.25 
Graphdiyne Present work a 167.26 6.99 2.37 0.21 1.79 
b 167.25 7.03 2.37 0.22 1.92 

Ref. 14-16 a 158.57 6.30 2.09 

b 144.90 6.11 2.19 


shown in Table 2. It is found that the relaxation times for 
electrons are longer than those for holes in both y-graphyne and 
graphdiyne. As is known, the deformation potential constant 
characterizes the strength of charge-phonon coupling, and 
weaker charge-phonon coupling leads to longer carrier relaxation 
times. The calculated Epp and t suggest that n-type doping may 
bring better thermoelectric performance than p-type doping in 
the two kinds of graphyne sheets. 

According to eqn (2)-(5), the electrical conductivity o, the 
electronic thermal conductivity Ke and the power factor S’o are 
calculated with respect to the electron relaxation time t. As 
shown in eqn (6), carrier relaxation times of two-dimensional 
sheets are inversely proportional to temperature. Within this 
relationship, the temperature-dependent relaxation times are 
determined based on the room temperature values, as summarized 
in Table 2. Fig. 3 plots the electronic transport coefficient map 
as a function of temperature T and carrier concentration n for 
y-graphyne (the upper panels) and graphdiyne (the lower 
panels). In Fig. 3 the positive carrier concentration corresponds to 
n-type doping while the negative one corresponds to p-type doping. 

It can be seen from Fig. 3(a) that the Seebeck coefficient S 
exhibits two obvious peaks near the Fermi level (n = 0), and the 
peaks of graphdiyne are higher than those of y-graphyne owing 
to the slightly larger band gap and doubly-degenerate bands. As 
to the electrical conductivity o shown in Fig. 3(b), there is a 
sharp increase with increasing doping concentration when 
temperature is considered as constant. The electrical conductivity 
o of y-graphyne is higher than that of graphdiyne. The obtained o 
for n-type doping is much larger than that for p-type doping, 
which is consistent with the above discussions about carrier 
relaxation times. As is known, higher o requires higher carrier 
concentration whereas higher S requires lower carrier concentration. 
Thus an optimized power factor S”o requires a certain compromise 
between S and a. As shown in Fig. 3(c), the maximum $’o appears 
at a moderate doping level where neither S nor ø reaches its 
maximum. The calculated electronic thermal conductivity 
(Fig. 3(d)) shows a similar behavior to that of electrical conductivity 
because of the Wiedemann-Franz law.** With higher Seebeck 
coefficient and lower electrical conductivity, graphdiyne exhibits 
higher power factor S*o and lower Ke than y-graphyne at same 
doping levels. 

If we focus on the high temperature areas (>450 K) of 
Fig. 3(c), we may find that the optimized power factor decreases 
with the increase of temperature when the carrier concentration 


(a) S(VK") (b) o(Sm") 


0 l 
C 600 
x 
E 500 
300 
-8 4 0 4 8 
a} 
-4 0 4 8 


n (10 cm*) 


800 


S 600 
E 500 
400 


300 


4 
n (10° cm?) 


0 4 8 


Fig. 3 


is considered as constant. In contrast, the electronic thermal 
conductivity slightly increases with the increase of temperature, as 
shown in Fig. 3(d). Thus the ZT values of these two graphyne sheets 
will reach their maximum near moderate temperature (400-900 K). 
In addition, compared with the power factor, the electronic thermal 
conductivity of graphyne increases faster with increasing carrier 
concentration, ranging 10'°-10*° cm7~*. Such a trend indicates 
that y-graphyne and graphdiyne could be optimized to exhibit 
the highest ZT values at a relatively low carrier concentration. 


3.3 Phonon dispersion relations and density of states 


Fig. 4 displays the calculated phonon dispersions for graphene, 
y-graphyne and graphdiyne. The gradually changed colors are 
determined by the group velocity, and the red and blue color 
stand for large and small group velocity, respectively. In general, 
our calculated phonon spectrum of graphene is consistent 
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Calculated electronic transport coefficient map as a function of temperature T and carrier concentration n for y-graphyne (top) and graphdiyne (bottom). 


with the previous experimental measurements and theoretical 
calculations.*°*° As shown in Fig. 4(a), the largest acoustic 
frequencies of graphene are almost as high as 1200 cm *, and 
the acoustic branches disperse with very high group velocity. 
Compared with graphene, two kinds of graphyne sheets exhibit 
much lower acoustic branches (<300 cm~') and introduce some 
new optical branches at around 2200 cm‘. The optical branches 
between 200 and 1500 cm * are fairly flat. Both y-graphyne and 
graphdiyne exhibit lower acoustic group velocity than graphene. 
The obtained phonon spectrum of y-graphyne is consistent 
with the previous DFT calculations.*”°** In the case of graph- 
diyne, the Raman spectra were recorded by Li et al. They found 
that the Raman peak at 1382.2 cm‘ is attributed to the 
breathing vibration of the sp? carbon in aromatic rings while 
the one at 2189.8 cm * is attributed to the vibration of the sp 
carbon in acetylenic linkages.® Our calculated phonon spectrum 
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Fig. 4 Calculated phonon dispersion relations for (a) graphene, (b) y-graphyne and (c) graphdiyne. The gradual change of colors from blue to red 
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Fig. 5 Calculated phonon density of states (DOS) for graphene, y- 
graphyne and graphdiyne (a) and sp-, sp-carbon partial density of states 
(PDOS) for (b) y-graphyne and (c) graphdiyne. 


of graphdiyne confirms the observed two phonon branches at 
1339.8 cm” * and 2216.5 cm’. 

Fig. 5(a) presents the phonon density of states (DOS) for gra- 
phene, y-graphyne and graphdiyne. In graphene, the low-frequency 
(<900 cm’) portion of the phonon DOS is mainly contributed by 
the acoustic phonon. In graphyne, the contribution of the acoustic 
phonon shrinks to a narrow frequency range (<300 cm“’). The 
phonon DOS of graphyne sheets presents several peaks at low 
frequency, showing the decreasing group velocity of the acoustic 
phonon. The intermediate-frequency portion of the phonon DOS is 
also split into some neighboring peaks, which is consistent with the 
flat optical branches. The isolated high-frequency DOS around 
2200 cm * is consistent with the isolated optical branches. 

As mentioned before, graphyne contains sp and sp” hybri- 
dized bonds while graphene is all-sp” bonded. The sp and sp” 
carbon partial densities of states (PDOSs) for y-graphyne and 
graphdiyne are separated and shown in Fig. 5(b) and (c). It can be 
seen that the isolated high-frequency phonons (2100-2300 cm *) 
are completely contributed by sp carbon atoms, which is very 
much consistent with the previous experimental conclusion." 
The low-frequency phonons (<600 cm *) are mainly contributed 
by sp carbon atoms while the intermediate-frequency branches 
(600-1500 cm™*) are mainly contributed by sp? carbon atoms. As 
is shown, the phonon DOS of sp*-hybridized aromatic rings in 
graphyne has almost the same shape as that in graphene. Along 
the series of graphene, y-graphyne, and graphdiyne, the sp” 
component in lower-frequency phonons is gradually replaced 
by the sp one. The phonon DOS shows that the presence of the 
sp bonding decreases the group velocity of acoustic branches and 
introduces some new optical branches, implying much lower 
phonon thermal conductivity in graphyne than in graphene. 


3.4 Phonon thermal conductivity 


As discussed in Introduction, the phonon thermal conductivity of 
graphene and graphyne has been investigated by using the NEGF 
method within the ballistic limit.*’? ?*”? The NEGF results show a 
monotonously increase of phonon thermal conductance with 
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Fig. 6 Calculated phonon thermal conductivity as a function of tempera- 
ture for graphene, y-graphyne and graphdiyne. The fitting curves of Kph ~ 
1/T* are also shown. 


increasing temperature. In the NEGF method, the phonon transport 
is treated as ballistic transport, and the phonon-phonon interaction 
is neglected. This approximation is applicable at just low tempera- 
ture when phonon scattering is weak. From the room temperature 
to intermediate temperature range, phonon-phonon scattering 
processes dominate the phonon transport and should be reasonably 
treated. In this work, the adopted phonon BTE method is applied to 
handle the phonon scattering effect and describe the relationship 
between phonon thermal conductivity and temperature. 

Fig. 6 displays the calculated phonon thermal conductivity as a 
function of temperature for graphene, y-graphyne and graphdiyne. 
The «pn of graphene at room temperature is 3845 Wm ' Kt, 
which is consistent with the experimental measurement (2500- 
5000 W m * K *”*) and close to the previous theoretical calcula- 
tions (3435-3800 W m * K *°*”). As shown in Fig. 6, our results are 
roughly scaled by x,,1/T* (@ = 1.66, 0.93, and 0.75 for graphene, 
y-graphyne and graphdiyne, respectively). The decreased thermal 
conductivity at high temperature can be owed to the more 
frequent phonon-phonon scattering. The calculated values for 
these two graphyne sheets are much lower than that for gra- 
phene in the whole temperature range. At room temperature, 
the thermal conductivities are 82.3 and 22.3 Wm‘ K* for 
y-graphyne and graphdiyne, respectively. These results are con- 
sistent with the above analysis of phonon dispersion and DOS. 

In a quantitative comparison, the phonon thermal conductivity 
of y-graphyne and graphdiyne is only 3% and 1% of that of their 
allotrope graphene, respectively. It is meaningful to understand 
the huge difference in the phonon transport properties between 
these similar two-dimensional materials. As is known, phonon 
thermal conductivity of two-dimensional sheets can be written as: 
Kph = Cybgl/2, where Cy, ¥g, and 7 are the specific heat, the group 
velocity, and the phonon mean free path (MFP). In this work, the 
phonon mean free path 7 is calculated as the the average of free 
path over all the phonon modes: 
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where Nqs is the number of g-points sampled. 


Table 3 Calculated specific heat Cy, mean values of group velocity Vg, 
phonon mean free path I, and thermal conductivity Kpn for graphene, y- 
graphyne and graphdiyne at 300 K 


Structure Cy (J cm™° K~) 9, (km s™*) 7 (nm) Kph (W m™ K7*) 
Graphene 1.64 5.42 866 (775*°) 3845 
y-Graphyne 1.68 1.65 59.5 82.3 
Graphdiyne 1.47 1.66 18.3 22.3 


The calculated Cy, v,, and 7 at 300 K for graphene, 
y-graphyne and graphdiyne are summarized in Table 3. As is 
shown, there is no obvious difference between the specific heat 
values while the mean group velocity of graphene is almost 
triple of that of graphyne. Our calculated phonon MFP for 
graphene is 866 nm, which is very close to the experimentally 
measured value of 775 nm.** The phonon MFPs of y-graphyne 
and graphdiyne are nearly one or two orders shorter than that 
of graphene. Considering the shorter phonon MFP, the much 
lower phonon thermal conductivity of graphyne can be under- 
stood. The previous NEGF calculations assumed that graphyne 
and graphene have the same MFP of 775 nm, and their 
obtained phonon thermal conductance op, of graphyne is 
60% lower than that of graphene.”*? Our studies show that 
the mixed sp/sp” bonds in graphyne sheets increase phonon- 
phonon scattering, and the calculated phonon MFP of graphyne 
is much shorter than that of graphene, resulting in the two 
orders lower thermal conductivity. 

It is interesting to examine the distribution of phonon free 
paths in graphene, y-graphyne and graphdiyne. Fig. 7 presents 
the normalized cumulative phonon thermal conductivity as a 
function of the free path at 300 K. Their phonon thermal 
conductivity distribution curves look very similar, except for 
the different scales of phonon free paths. Usually, low-frequency 
phonon modes with large free paths dominate in the anharmonic 
scattering processes, and high-frequency modes with small free 
paths contribute a little to the thermal conductivity.°”°° As shown 
in Fig. 7, the thermal conductivity of graphene is mostly contributed 
from long free path phonons, and the fraction of heat transport 
carried by phonons with free paths longer or equal to 1000 nm is 
~55%. The phonon free paths of y-graphyne are distributed in 


the range of 1-5000 nm, and the longer free path phonons 
(>100 nm) nearly contribute half of the thermal conductivity. 
The free paths of graphdiyne are almost shorter than 1000 nm. 
Its thermal conductivity is mostly contributed by short free path 
phonons, and the fraction of heat carried by a phonon with the 
free path shorter or equal to 20 nm is ~55%. 


3.5 Optimized ZT values 


With the obtained transport coefficients, we are now able to 
evaluate the thermoelectric performance of graphyne according 
to ZT = S°oT/(Kk. + Kpn). Fig. 8(a) and (b) presents the ZT value as 
a function of temperature T and carrier concentration n for 
y-graphyne and graphdiyne, respectively. The red and blue 
colors stand for large and small ZT values, respectively, as 
shown in the color scale. As can be seen from Fig. 8, the calculated 
ZT values exhibit highlighted areas in a wide temperature range. 
This observation suggests that one can always obtain high thermo- 
electric performance of graphyne by appropriate doping. The 
n-type graphyne sheets exhibit larger ZT values than the p-type 
systems, which is mainly caused by the relatively high electrical 
conductivity and the power factor of n-type doping (see Fig. 3(b) 
and (c)). In addition, due to the higher power factor (see 
Fig. 3(c)), lower electronic thermal conductivity (see Fig. 3(d)) 
and lower phonon thermal conductivity (see Fig. 6), graphdiyne 
can be doped to exhibit larger ZT values than y-graphyne. 

The optimized ZT values of y-graphyne and graphdiyne at 
300 K are summarized in Table 4, where the corresponding 
carrier concentration n and transport coefficients S, 0, S?o, Ke, 
and «pn are also given. Wang et al.° reported that the optimized 
ZT of graphene at room temperature is 0.0094. The corresponding 
transport coefficients are extracted from their studies and listed in 
Table 4 for comparison. It is found that graphyne could exhibit a 
higher power factor than graphene by light doping (~10'? cm~?). 
As shown in the previous reports, the mixed sp/sp* bonds open a 
band gap in graphyne sheets.‘ ”* The low electronic thermal 
conductivity, the high Seebeck coefficient and the improved power 
factor in graphyne can be owing to the presence of sp bonding. On 
the other hand, sp bonding in graphyne sheets leads to low 
thermal conductivity. Reasonably, high ZT values are obtained 
for graphyne sheets. Compared with y-graphyne, graphdiyne 
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Fig. 7. Normalized cumulative phonon thermal conductivity as a function of the mean free path at 300 K for (a) graphene, (b) y-graphyne and (c) 


graphdiyne. 
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Fig. 8 Calculated ZT value map as a function of temperature T and carrier concentration n for (a) y-graphyne and (b) graphdiyne. 


Table 4 Optimized ZT values of y-graphyne and graphdiyne at 300 K as well as the optimal operating temperature. The corresponding optimal carrier 
concentration n, the Seebeck coefficient S, the electrical conductivity a, the power factor S?s, and the electronic and phonon thermal conductivity xe and 


Kpn are also indicated 


Structure n (10*° cm ~*) S (uV K+) o (10° S m~’) 
Graphene? @300 K ~77 ~16.4 
y-Graphyne —4.0@300 K 260 5.55 
2.1@300 K —316 8.96 
—5.5@780 K 286 2.78 
3.3@760 K —307 5.76 
Graphdiyne —4.8@300 K 248 1.91 
1.1@300 K —372 4.09 
—7.0@800 K 315 0.95 
1.2@580 K —429 2.34 


contains more sp bonding, and it exhibits lower thermal con- 
ductivity and higher ZT value. 

As shown in Table 4, our calculations indicate that the two 
orders lower phonon thermal conductivity in graphyne sheets 
significantly increases the ZT values. The optimized ZT value of 
graphdiyne at room temperature is 0.99 and 3.26 by p-type and 
n-type doping, respectively. These values mostly exceed the 
previous reports and make graphyne a very promising candi- 
date for thermoelectric applications. In the previous studies, 
Wang et al. and Sevinçli et al. independently studied the 
thermoelectric properties of y-graphyne sheets by using the 
NEGF method within the ballistic approach. They obtained 
similar results and the optimized ZT value for y-graphyne is 
only about 0.16 at 300 K.””° They attributed the relatively large 
ZT values of graphyne to larger Seebeck coefficients. The huge 
difference in the phonon mean free path between graphene and 
graphyne is neglected by the NEGF method, which leads to 
overestimated thermal conductivity and underestimated ZT 
value for graphyne sheets.””* Sun et al. investigated the thermo- 
electric properties of graphdiyne by a combination of the BTE 
method and MD simulations, and the optimized ZT value for 
graphdiyne is 4.8 at 300 K.” As is known, MD results strongly 
depend on the empirical potential coefficients. It seems that the 
thermal conductivity of graphene and graphyne obtained by the MD 
method is underestimated compared with experimental values. For 
example, the MD results show that the lattice thermal conductivity of 
graphene at 300 K is only 100 W m7? K7','8 which is much smaller 


than the experimental measurements (2500-5000 W m~t K_*).°” 


S’o (Wm! K~?) Ke (W m™* K~’) Kph (W m+ K?) ZT 
~0.097 ~130 ~ 2960 0.0094 
0.38 40.8 82.3 0.91 
0.89 65.9 82.3 1.81 
0.23 53.0 34.4 2.02 
0.55 107.2 35.2 2.92 
0.12 14.0 22.3 0.99 
0.57 30.1 22.3 3.26 
0.09 18.7 10.8 2.56 
0.43 33.2 13.8 5.30 


The high ZT value of 4.8 for graphdiyne reported by Sun 
may result from the underestimated thermal conductivity 
(7.3 Wm * K “at 300 K). 

It can be seen from Fig. 8 that y-graphyne and graphdiyne 
exhibit very high ZT values in a wide operating temperature 
range around 800 K and 600 K, respectively. The highest 
ZT values and the corresponding transport coefficients are 
summarized in Table 4. Compared with the results at room 
temperature, the optimized ZT values increase from 1.81 and 
3.26 to 2.92 and 5.30 for n-type y-graphyne and graphdiyne, 
respectively. The increased ZT values result from the exponential 
decline of phonon thermal conductivity and the mild decrease of 
the power factor. As shown in Table 4, the electronic thermal 
conductivity is comparable and even in excess of the phonon 
contribution, which is not uncommon in doped systems.’ We 
note that the ZT value of graphdiyne is higher than 3.0 in a very 
wide temperature range (almost 300-900 K) and a wide doping 
range (n = 0.2-8.0 x 10‘? cm-*), which implies that this 
graphdiyne is particularly suitable for fabrication of the thermo- 
electric modules. 


4 Summary 


In summary, the electronic, thermal, and thermoelectric trans- 
port properties of y-graphyne and graphdiyne are investigated 
by using first-principles calculations and the Boltzmann trans- 
port equation method. Graphyne exhibits much higher Seebeck 


coefficients than graphene, which is consistent with its semi- 
conducting band structure. The deformation theory is adopted 
to describe the electron-phonon interaction and our results 
predict that the carrier relaxation times of electrons are larger 
than those of holes, implying better thermoelectric perfor- 
mance of n-type graphyne. Contrary to the previous NEGF 
calculations, our calculated lattice thermal conductivity of 
graphyne is found to be two orders lower than graphene. The 
phonon spectra and DOSs indicate that the low lattice thermal 
conductivity originates from the decrease of the group velocity 
and the appearance of new phonon branches. The calculated 
phonon MFPs for graphene, y-graphyne, and graphdiyne are 
866, 60, and 18 nm, respectively. We argue that the overestimation 
of lattice thermal conductivity in previous studies is due to the 
rough description of phonon-phonon scattering. With the high 
power factor and relatively low thermal conductivity, graphdiyne 
possesses superior thermoelectric performance with a maximum 
ZT value of 5.3 at 580 K. Our study suggests that major effort 
should be paid to the realization of these high-performance and 
environmentally-friendly thermoelectric materials. 
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